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We use soft collinear effective field theory (SCET) to study nonperturbative strong interaction 
effects in Z decays to hadronic final states that are enhanced in corners of phase space. These 
occur, for example, in the jet energy distribution for two jet events near Ej — Mz/2, the thrust 
distribution near unity and the jet invariant mass distribution near zero. The extent to which such 
nonperturbative effects for different observables are related is discussed. 
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I. INTRODUCTION 

Some of the most successful applications of perturba- 
tive QCD are to processes such as Z decay to hadrons 
or e + e~ annihilation at large center-of-mass energy, in 
which a state with no strong interactions decays into final 
hadronic states. This paper will discuss the case of Z de- 
cay, but the results in this paper apply equally well to the 
other cases. Not only is the total hadronic Z decay width 
calculable but so are less inclusive infrared-safe quantities 
like the Z decay rate into 2-jet and 3-jet events, the thrust 
distribution and jet mass distributions. Comparison of 
perturbative predictions for these and other quantities 
with experimental data on Z decays from LEP and SLD 
has led to a remarkably accurate extraction of the strong 
coupling constant a s {M z ) 0, H S H H H . Although 
the extraction of the strong coupling from event shape 
variables is less accurate than from the total hadronic 
Z width, it is more model-independent since (neglecting 
quark mass effects) it does not depend on the values of 
the quark couplings to the Z. 

For the totally inclusive hadronic Z decay width, the 
operator product expansion allows one to include in theo- 
retical predictions nonperturbative strong interaction ef- 
fects that are characterized by vacuum expectation values 
of local operators. The effects of higher-dimension oper- 
ators are suppressed by powers of the strong interaction 
scale Aqcd divided by the center-of-mass energy Mz- 
Since the Z mass is large, these effects are very small. 
For example, if quark masses are neglected, the leading 
nonperturbative effects in the Z decay width come from 
the vacuum expectation value of the gluon field strength 
tensor squared, (G /J „G A " / ). This dimension-four opera- 
tor gives rise to corrections to the total hadronic width 
suppressed by Aq CD /M| ~ 10~ 9 . 

Less inclusive variables that characterize Z decay to 
hadrons give rise to nonperturbative effects suppressed 
by smaller powers of A QC d/M z [1 @, 111 13 
Ibi . Il5j . Furthermore, these corrections often become 
even more important in corners of phase space where 
hadronization effects are significant, such as in the thrust 
distribution very near T = 1. It has been conjectured 
that the enhanced nonperturbative effects to many event 
shape distributions have a universal form with a single 
nonperturbative parameter 0,0,0,0,0,0. These 



arguments are based on analysis of renormalon ambigu- 
ities in the QCD perturbation series and on the behav- 
ior of rcsummed perturbation theory. The conjectured 
relationship between the nonperturbative corrections to 
event shape distributions has recently been tested exper- 
imentally |2j. 

Recently the enhanced nonperturbative effects that oc- 
cur for the jet energy distribution in corners of phase 
space have been studied using effective field theory meth- 
ods 0| . This approach uses the fact that very low mo- 
mentum degrees of freedom which contain the nonper- 
turbative physics couple to the degrees of freedom with 
energies of order Mz via Wilson lines. Nonperturbative 
effects have been extensively studied previously us- 
ing factorization methods to divide the process into hard, 
jet-like and soft subprocesses 0,0]. Nonperturbative 
effects are computed from the soft subprocess. The ef- 
fective field theory approach is similar to the one based 
on factorization methods. In this paper we elaborate on 
the work in [19j and extend it to other shape variables. 
The enhanced nonperturbative effects are expressed in 
terms of weighted matrix elements of operators involving 
Wilson lines, where the weighting depends on the event 
variable being considered. We hope that this paper will 
help make the results of Ref. 0] more accessible to the 
community of high energy theorists who are most familiar 
with effective field theory methods. 

In this paper we study smeared distributions which al- 
lows us to expand the nonperturbative effects in powers 
of Aqcd, and write them as matrix elements of Wilson 
line operators and their derivatives. The computations 
are similar to those of smeared distributions in the end- 
point region in B decay — the point-by-point computa- 
tion requires knowing the nonperturbative shape func- 
tion, whereas nonperturbative effects in the smeared dis- 
tributions can be written in terms of Ai ; 2 provided the 
smearing region is large enough. 

For pedagogical reasons we start with a detailed treat- 
ment of the jet energy Ej in Z decay to two jets, where 
the jets are defined as Sterman and Weinberg did in their 
original work on jets in QCD 0] . We spend considerable 
effort on this variable because the theoretical expression 
for its enhanced nonperturbative corrections is simpler 
than for other more phenomenologically interesting vari- 
ables like thrust. Cone algorithms for jets, like that of 
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Sterman and Weinber g ar e ambiguous at higher orders 
in perturbation theory |23l |24| . However, since our main 
interest is in nonperturbative effects and we only treat 
the perturbative effects at order a s , that ambiguity is 
not important for us. At lowest order in perturbation 
theory, the Z boson creates a quark and an antiquark, 
each with energy Mz/2, and so the jet energy distribu- 
tion is equal to 



dr 



2-jct 



dEj 



= r 



(0) 
2-jct 



6(Ej - Mz/2), 



(1) 



where r^j et is the total two-jet rate at lowest order in 
perturbation theory. This leading order theoretical ex- 
pression for the jet energy distribution is singular at 
Ej = Mz/2. Furthermore, the leading perturbative 
and nonperturbative corrections are also singular at that 
kinematic point. However, a non-singular quantity that 
can be compared with experiment without any resum- 
mation of singular terms is obtained by smearing the jet 
energy distribution over a region of size A that contains 
the lowest order partonic endpoint at Ej = Mz/2. The 
leading nonperturbative correction to this smeared en- 
ergy distribution is suppressed by Aqcd/A. So, for ex- 
ample, with A ~ 10 GeV the nonperturbative corrections 
are expected to be of order 10%, roughly the same size 
as perturbative corrections, and an order of magnitude 
larger than the order Aqcd/^z correction expected in 
the complete two jet rate. We argue that for Ej very near 
Mz/2 it is not possible to capture the dominant nonper- 
turbative effects simply by shifting, Ej — > Ej — /i np , in 
the perturbative expression for dT2-jet/<iEj (where /x np 
is a nonperturbative parameter of order Aqcd)- 

In the next section, we derive an expression for 
the leading enhanced nonperturbative correction to the 
smeared jet energy distribution for two jet events us- 
ing methods from soft-collinear effective field theory 
(SCET) EHH. This correction is given by the 

vacuum expectation value of a nonlocal operator involv- 
ing Wilson lines. Perturbative order a s corrections to 
this variable are derived in Appendix iBl 

Section Infl discusses the leading nonperturbative cor- 
rections for thrust, jet masses, the jet broadening vari- 
ables, the C parameter and energy-energy correlations. 
In agreement with Ref. 0] we nn d that the correction to 
jet mass sum and thrust are related. However, without 
additional model-dependent assumptions we do not find 
that the enhanced nonperturbative corrections to the C 
parameter and jet broadening variables can be related 
to those for thrust and the jet masses. We compare the 
level of our understanding of the enhanced nonperturba- 
tive effects in these variables. 



II. OPERATOR PRODUCT EXPANSION FOR 
THE TWO JET ENERGY DISTRIBUTION 

The nonperturbative corrections to the energy distri- 
bution for Z decay to two jets, dr2_j e t/d-Ej near Ej = 



Mz/2 are computed in this section. The perturbative 
corrections will be discussed in Appendix IBl The results 
are given for the Sterman- Weinberg jet definition, where 
a cone of half-angle S contains a jet if the energy con- 
tained in the cone is more than E cut — (3Mz- We take 
the cone half-angle 5 and the dimensionless energy cut 
variable (3 to be of order a small parameter A, and com- 
pute in a systematic expansion in powers of A. We are 
interested in the jet energy distribution within a region 
A of Mz/2, where M z > A > X 2 M Z . For example, 
A ~ XMz- The Sterman- Weinberg jet definition, like 
other cone algorithms, is ambiguous. However, this diffi- 
culty does not become apparent at order a s in perturba- 
tion theory and also will not influence our discussion of 
nonperturbative effects. 

SCET is the appropriate effective field theory for 
the kinematic region of interest, and will be used for 
the derivation of the nonperturbative corrections to 
dr2_j c t/d£',/ near Ej = Mz/2. It is convenient to 
introduce two lightlike vectors n and n which satisfy 
n° = n° = 1 and n = — n. Four-vectors are de- 
composed along the n, n and perpendicular directions: 

V = (V + ,V~,Vj_) where V + — n ■ V, V~ — n ■ V and 

VI = V- V+n^ /2 - V-n» /2. For the problem of inter- 
est, SCET contains n-collinear, n-collinear and ultrasoft 
degrees of freedom j^. The n-collinear and n-collinear 
degrees of freedom have typical momenta that scale as 



(") ~ M Z {\\ 1, A), pf> ~ M z (l, A 2 , A), (2) 



P 



and the ultrasoft degrees of freedom have momenta that 
scale as 



p u ~M z (\ 2 ,\ 2 ,\ 2 ). 



(3) 



We take A ~ \/ Aqcb/Mz which implies that the typ- 
ical "off-shellness" of the ultrasoft degrees of freedom, 
Pu ~ M§ A 4 ~ Aq CD , is set by the QCD scale while 
the typical "off-shellness" of the collinear degrees of free- 
dom, p 2 ~ M| A 2 ~ Mz Aqcd , is much larger than Aq CD . 
Hence the collinear degrees of freedom can be treated in 
perturbation theory. 

The nonperturbative effects we are after are charac- 
terized by matrix elements of operators composed from 
the ultrasoft degrees of freedom. In Z decay into two 
jets, the jets are almost back-to-back, and n is chosen 
along one of the jet directions. The degrees of freedom 
in the two jets are then represented by n-collinear (for the 
antiquark jet) and n-collinear fields (for the quark jet). 
In this section we work to lowest order in perturbation 
theory in the collinear fields. Hence we match the weak 
neutral current in full QCD onto the effective theory at 
tree level, 



3 



(4) 



where — gv7± + 9Al±"/5 involves the vector and ax- 
ial couplings of the Z boson. The fields £ s and are 
collinear quark fields in the n and n directions and we 
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have adopted the convention 



p 



(5) 



where the label momentum p contains the components 
of order 1 and A, n ■ p and p±, and the order A 2 compo- 
nents are associated with the space-time dependence of 
the fields. The Wilson lines W n ,n are required to ensure 
collinear gauge invariance |27| . Since in this section we 
work to lowest order in QCD perturbation theory, they 
play no role in the analysis and can be set to unity. 

The typical momenta of the partons in the jets are of 
the order of the collinear momenta, Eq. where the 
overall scale of their momentum is set by Mz- However, 
it is possible for the jets to contain partons with momenta 
that have an overall scale that is much less than Mz- 
Because of the sum over all values of p in Eq. Q, such 
partons can still be represented by collinear fields. The 
interaction of n-collinear fields among themselves is given 
by the full QCD Lagrangian, and so the hadronization 
of n-collinear partons into a jet is described by the full 
theory. 

The Lagrangian of the effective theory does not contain 
any direct couplings between collinear particles moving in 



the two different lightlike directions labeled by n and n; 
however, they can interact via the exchange of ultrasoft 
gluons. It is convenient to remove the couplings of the 
collinear degrees of freedom to the ultrasoft ones via the 
field redefinition [H 



Cn y ^n^-rt: A n ► Y^J A n Y n , 



where A n is an n-collinear gluon field and 



Y n (z) = Pexp 



ig j ds n ■ A u (ns + z) 
o 



(6) 



(7) 



denotes a path-ordered Wilson line of ultrasoft gluons in 
the n direction from s — to s = oo. This is the appro- 
priate field redefinition for outgoing collinear fields, since 
if a factor of exp(— es) is inserted in the integrand to de- 
couple the interactions at late times, one reproduces the 
correct it prescription for the collinear quark propaga- 
tor. For annihilation which contains incoming collinear 
particles Y n is from s = — oo to s = and the daggers 
are reversed in Eq. ©. An analogous field redefinition 
with n — y n removes the couplings in the Lagrangian of 
ultrasoft fields to the n-collinear fields. 



The differential decay rate for Z decay to two jets is 
1 

2M~z ., 3 



dr 



2-jct 



E \(JnJnX u \f(0)^ |0)| 2 {2*)H\p z - PJn ~p Jn - ku) 



final states 



where the sum over final states includes the usual phase space integrations and e is the polarization vector of the 
decaying Z boson. Since after the field redefinitions shown in Eq. ffijl. there are no interactions between the ultrasoft 
and collinear degrees of freedom, the matrix element factorizes, and at lowest order in perturbation theory in the 
collinear degrees of freedom, 



dr 



i 



2-jct 



d 3 P 9 



d 3 P 5 



2M Z (27r) 3 2p0 (27r) 3 2p? 
1 



M 



(o) 

if 



^2{2ir) 4 6 4 (p z -p q -Pq - k u ) 



T[Y nd e Yl a }(0) X u (k u )) (X u (k u ) T[Y na C Y^ C d ](0) 



(9) 



In Eq. (O, \-M-\f | 2 is the square of the Z — » qq decay am- 
plitude averaged over Z polarizations and summed over 
the quark and antiquark spins and colors, T(T) denotes 
time (anti-time) ordering, Nq is the number of colors, 
and we have explicitly displayed the color indices on the 
ultrasoft Wilson lines. 

The derivation of Eq. @ in many ways parallels the 
use of the operator product expansion to compute the 
deep inelastic scattering cross-section, or the rate for in- 
clusive semileptonic B decay. There is, however, one im- 
portant distinction. The sum over final states in deep 
inelastic scattering and B decay is a sum over a com- 
plete set of color singlet hadron states. In Eq. (jSJ), one is 



summing over a complete set of jet and ultrasoft states. 
These are a complete set of partonic states, and are not 
necessarily color singlet states. In fact, unitarity would 
be violated if one separately imposed the color singlet 
condition on each of \J n ), \Jn) and \X U ). The deriva- 
tion of Eq. is valid to the extent that the sums over 
partonic and hadronic states are equivalent. In jet pro- 
duction, the color of the fast quark that turns into a jet 
is eventually transferred to low energy partons during 
the fragmentation process. The low energy partons com- 
municate between the different jets, and make sure the 
whole process is color singlet. The assumption is that 
this color recombination does not affect the decay rate 
at order Aqcd/Mz. 
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To calculate dlYjet/d-Ej we integrate Eq. © over 
the allowed values of the quark and antiquark three- 
momentum with the factor 5{Ej — p q ) inserted. This 
corresponds to choosing the quark jet as the "observed" 
jet. If one does not distinguish between quark and an- 
tiquark jets then Eq. still applies since the value of 
d^-jet/d-Ej when the "observed" jet is an antiquark jet 
is the same. It is convenient to work in the rest frame of 
the decaying Z, pz — (Mz,Mz,0±), and align n with 
the quark three-momentum p q . The decomposition of 
the quark's four-momentum in terms of label and residual 
momentum, p q = p q +k q , has the formp+ = p^ + k^ with 
Pq = 0, p q ± = 0. (Note this means that p q ± = k q ± = 
and k~ — 0.) Hence the phase space integration over 
quark three-momentum becomes 

Pi 

For the antiquark's four-momentum the decomposition 
into residual and label momentum is p^ — k^ , p^ — 



p^ + k^ and p q ± = p q ± + kq± . One cannot set p^_L = 
by a choice of n, since n = n, and n has already been 
fixed by the direction of the quark jet. 

Expressed in terms of label and residual momenta the 
phase space integration over antiquark three-momentum 
is 

Pq 

- 2^J 2(2tt)3 p • 1 ' 

Pq 

Here the delta function S((pq + kg) 2 ) = 5{pqk^ — p|j_) 

was used to do the kq integration setting k^ = Pq±/pq ■ 
At leading order in the SCET expansion parameter A 
the invariant matrix element M-if only depends on the 

label momenta p+ and pq . In terms of label and resid- 
ual momentum the energy-momentum conserving delta 
function becomes 



S 4 (pz - p q - Pq - k u ) = 2S(p z -pq - k u )8{p% -P+ -p% - k+)8 2 {p q ±_ +k n _L) 

= 2S M z ,p- S M z ,p+ S l^,oHkq + K)S(k+ + K)S 2 {^± + K±)- 



(12) 



The relation kq = p 2 q ^_/pq and the Kronecker delta that sets Pg_L to zero imply that kq = 0, and so this variable 
does not appear in the penultimate delta function in Eq. \12\ . 
Using these results gives, 



cir 



2-jet 



\M^ 2 



dE, 



8M Z (2 



L J d fe+ J dkgd 2 k ql _ J2s(^-E J+ ^pj 5{kz + k-)S(k+ + k+)S 2 (k ql _ + Kx) (13) 



(o| f[Y nd e Yl a ](0) \x u (k u )) (x u {k u )\ T[Y na c Yl d ](0) |o) 



\M 



(0)|2 
if 



WirMz 



7zJ2 s (^f- E J-T)w(° f[Ynd eyL a](0) Xu{ku) ) ( Xu{ku) T[Y ™ cyl d](0) 



We write the remaining delta function as the integral 



Mz „ K 



du 
2^ 



exp 



Mz „ K 



(14) 



At this stage the collinear degrees of freedom have been 
integrated out, and the matrix elements above, which 
involve only ultrasoft degrees of freedom, are evaluated 
at leading order in the SCET expansion parameter (i.e. 
A — ► Q). Recall that the Sterman- Weinberg jet crite- 
ria restricts particles outside the cones used to define 
the two jets associated with the quark and antiquark to 
have energy less than E cut which we are taking to be or- 
der XMz- In the limit A — > this energy cut becomes 



much larger than a typical component of an ultrasoft 
four-momentum. Hence, for the matrix elements of these 
operators, E cut should be taken to infinity and does not 
restrict these matrix elements. Similarly the cone angle 
is taken to be of order A while the typical angle between 
components of ultrasoft momenta is order unity. Thus 
the cone angle should be taken to zero in the effective 
theory that contains only ultrasoft degrees of freedom 
and so there is no restriction on the ultrasoft states that 
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are summed over in Eq. (|13(l . 

Using the exponential dependence on k u to translate 
the anti-time ordered product to the space-time point 
un/2, and then using the completeness relation to per- 
form the sum over all ultrasoft intermediate states, we 
find for the jet energy distribution 



dIVjct = if ] ct S(M z /2 - Ej) 



dE. 



j 



(15) 



where the shape function S is defined by ^(J 



(16) 



x (0 



T [Y nd e Yl a ]{un/2)T[Y na C Y^ C d }(0) 



and the total two jet Z-decay width at lowest order in 
perturbation theory is 



.(0) 
2-jet 



16irM z 



NgMz , 2 , 2 , 

— 77. — [9v +9a)> 
LZtt 



(17) 



having implicitly summed over spins and colors. The n- 
directed and n-dircctcd ultrasoft Wilson lines commute 



since (sin — s 2 «) = — 4siS2 < 0, and the gauge fields in 
the Wilson lines are space-like separated. 

In this derivation we chose the jets to be composed 
entirely of collinear degrees of freedom. This is appropri- 
ate since jets are confined to narrow cones. For example 
the momentum of any massless particle in the quark jet 
satisfies p~ <C p + , which is the appropriate scaling for 
collinear particles in the n direction. However, it is pos- 
sible to repeat the above derivation allowing ultrasoft de- 
grees of freedom to be inside a jet. Then instead of insert- 
ing S(Ej — p°) into Eq. Q, one inserts S(Ej — p® - k°j), 

where k°j = (fc+ 7 /2)[l + 0(A)] denotes the total ultrasoft 
energy inside the quark jet. Using the delta functions in 
Eq. i|12fl we obtain again Eq. (|13fl . with fc+ in the final 
delta function now denoting the total ultrasoft momen- 
tum outside the quark jet. However, as mentioned pre- 
viously, at leading order in A the cone angle of the jet 
shrinks to zero, and one recovers the previous result. 

It is possible to remove the time and anti-time order- 
ing completely in the definition of the shape function S. 
Using the results from Appendix IA II our expression for 
the shape function becomes 



S(k) = 



1 



du 



i k u 



N c J 2tt 
1 



Nc 



[Yl d Yl a \(unl<2)[Y fia ^Y n c ](0) 





[y[ d Yl a ]5(k + m -d/2)[Yaa C Y n 



(18) 



where the ovcrlinc denotes an anti-triplet Wilson line. 

Since in the kinematic region of interest Mz/2 — Ej is much larger than n ■ d acting on ultrasoft gauge fields it is 
appropriate to expand the delta function above which gives 

S{M z /2 - Ej) = S(M z /2 - Ej) + S'(M z /2 - Ej) <0| O x |0) + ^S"(M z /2 - Ej) (0| 2 |0) + . . . (19) 



where 



0„ 



1 



Y jYJ 



in ■ d 



Y- C Y 

I 



N c 



-Tr 



in ■ D 



(20) 



The simple form for the operators O m arises because the 
variable Ej is totally inclusive on the "unobserved" an- 
tiquark jet. 

The formula for dIVjct /dEj is 



dr 2-jct _ p (0) 
dEj ~ 2 "j ct 



8{M z /2 - Ej) 
+5'{M z /2 - Ej) (0| Oi |0) 



(21) 



The delta function term in Eq. (|19[) simply reproduces 
the leading perturbative formula for dT^ ct /dEj while 



the higher-order terms contain the effects of nonpertur- 
bative physics. The derivation presented here assumes 
the observed jet is the quark jet. A similar derivation in 
the case where the antiquark jet is observed gives opera- 
tors 



O u 



1 

Nc 



-Tr 



Y, 



D 



Y, 



(22) 



Since the vacuum expectation values of O m and O m are 
equal by charge conjugation, our results also hold in the 
case where one does not distinguish between quark and 
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antiquark jets. 

We define the matrix elements using dimensional reg- 
ularization with MS subtraction so that in perturbation 
theory the vacuum expectation values (0| O m |0) are zero. 

Note that O2 is a very different operator than 0\ so it 
is not possible to capture the effects of nonperturbative 
physics for \Ej — M z /2\ ~ \ 2 M Z 1 simply by taking the 
lowest order perturbative formula in Eq. i|21|) and shifting 
Ej by a nonperturbative parameter /x np , i.e. Ej -> Ej - 
/i np . This ansatz results in the shape function 

S(M z /2 - Ej) = 6(M z /2 - E,,) + S'(M z /2 - Ej)fi np 
+ \5"(M z /2 - Ej)& v + ... (23) 

where the series of derivatives of delta functions has 
coefficients that are simply related by (0| O m 0) = 
(0| 0\ |0) m , which is not correct. 

For \Ej - M z /2\ ~ X 2 M Z all terms in the series 
of Eq. (|f 9fl are equally important. However for \Ej — 
M z /2\ ~ A 3> \ 2 M Z the vacuum expectation value 
of 0\ provides the leading order Aqcd /A nonpertur- 
bative correction. In this kinematic region the shift 
Ej — » Ej — /x np , with /i np = (0|Oi|0), correctly captures 
the most important effects of nonperturbative physics. 

We have focused on nonperturbative effects that are 
enhanced in the region near Ej = M z /2. If one considers 
a variable like the average value of the jet energy over 
the entire allowed phase space, then there are sources of 
nonperturbative corrections that we have not considered. 

Using the results of Appendix IA 21 the operator 0\ 
in Eq. I|2(J|I can be expressed in terms of the gluon field 
strength tensor [lflj : 



III. ENHANCED NONPERTURBATIVE 
CORRECTIONS TO EVENT SHAPE VARIABLES 



There are a number of event shape distributions that 
are commonly studied in the literature. Conventionally, 
one defines a general event shape distribution der/de, 
where e is an event shape variable defined such that the 
region e — > corresponds to the two jet limit. Examples 
are e = I — T for thrust, e = B for jet broadening and 
e = C for the C parameter. Any event shape distribution 
in Z decay contains both perturbative and nonperturba- 
tive contributions. The perturbative effects can be com- 
puted as a perturbation series in a s (M z ). At leading or- 
der, only two-jet (i.e. qq) events contribute. Events with 
more hard partons are suppressed by powers of a s (M z ). 
In general, nonperturbative effects are suppressed by 
powers of Aqcd /M z , but in corners of phase space where 
e < 1 these nonperturbative effects become enhanced. 
Here we consider the region Aqcd <C M z e <C M z and 
focus on the enhanced nonperturbative contribution sup- 
pressed only by a single power of AQcn/(M z e). 

Working to leading order in a s (M z ), the dominant 
nonperturbative effects are corrections to the two-jet dis- 
tribution. Nonperturbative corrections to higher order 
processes are suppressed by additional powers of a s (M z ) . 
We will compute the enhanced nonperturbative correc- 
tion to some commonly measured event shape distribu- 
tions, just as we did for the jet energy distribution in 
Sec. |TI] Recall for the jet energy distribution the domi- 
nant nonperturbative correction came from expanding 



M, 



M z 
2 



E., 



(25) 



+5> -Ej)^ 



2 



Oi = ^[Yi (in ■ D)Y n ] 



(24) 



i T r 
2 



ig l ds Y^{z]s 1 Q)n^n v G llv Y ri {z\s 1 0) 




Oi in Eq. (|24|l vanishes if the ultrasoft gauge field is 
treated as a classical degree of freedom. Then the Wil- 
son lines in Eq. (|24|) are unitary matrices and the trace 
vanishes since the gluon field strength tensor is in the ad- 
joint representation. Note that the vacuum expectation 
value of 0\ can still be nonzero because of quantum ef- 
fects. Usually operators involving products of gluon fields 
require renormalization, however, it is straightforward to 
show that 0\ is not renormalized at one loop. 



1 More correctly the differential cross section aT2-jet/d-Ej smeared 
over a region A of energy (that contains Ej = Mz/2) with A of 
order \ 2 M Z . 



in Eq. II Kill to linear order in fc+ . The delta function from 
Eq. I|12|) sets /c+ = — fc+, and we therefore find 



dr (o) 

Q1 2-jct _ r (o) 
dEj 2 "j ct 

where 

<^ + > = E 



(M z 



sr-^-Ej\s / r-^-E 



M z 



(K) 



(26) 



Nr 



T 



l nd 1 ne 



(0) X u (k u ) 



x (XM\T [Y na c Yl d ] (0) |o) k+. (27) 

The jet energy distribution has the nice property that 
one can write as the vacuum expectation value of an 
operator involving Wilson lines of ultrasoft gauge fields 
[namely, Eq. I|24|) ]. For some shape variables this is not 
possible. However, expressions analogous to Eqs. (|26H27|) 
can be derived. 
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A. Thrust 

First we consider the thrust distribution dT/dT where 
the thrust T is defined by 



M z T 



max 
t 



El*-Pil 



(28) 



where t is a unit vector that defines the thrust axis. 
The maximum is taken over all possible directions of t, 
and the sum is over all final state particles. To the or- 
der we are working the thrust axis t can be set equal 
to the spatial part of the lightlike four-vector n used 
to define the collinear antiquark field. It is convenient 
to call this direction the z-axis. The thrust distribu- 
tion is calculated analogously to the two jet distribu- 
tion except that the delta function 8(Ej — p®) is re- 
placed by S (M Z T -\Pg\- \p z q\ - J2 a \Ka\)i where the 
sum is over all ultrasoft particles. We adopt the same 
conventions as in the jet energy distribution so that the 
phase space integrals are again done using the delta func- 
tion in Eq. (|12fl . Decomposing the total ultrasoft four- 
momentum, k u = + ku^ , into the sum of the ultrasoft 
momentum from particles in the same hemisphere as the 
antiquark (type a) and the same hemisphere as the quark 
(type b) the thrust T can be written as 



M Z T 



Pa + Pa 



OL 

1 , 



= \r q +\k+\ (k - 4) + \k 

+1 (M- - k^+\ - - (k^~ - k^+\ 
+ 2 v / 2 V / ' 

Now the delta functions in Eq. I|12l) set p^ = p+ = Mz, 

kq = —k~, k+ = —k+, and kq — 0. Thus we find 



1 - 



M z \ u + u ) 



(30) 



where we have also used k u = ku + ku^ ■ Thus, 



dT 
dT 



,(o) 

2-jet 



.(0) 
2-jet 



8(1 — T) — 8'(1 — T) 



.(«)+ 



M, 



8(1-T)-S'(1-T) 



(Of) 



Af, 



(31) 



The thrust axis and the hemispheres are determined by 
the jet directions, and can be defined in terms of the la- 
bel momenta of the quark and antiquark. Thus t and 
the hemispheres a and b are label variables. Neverthe- 
less, because of the hemisphere condition on the ultrasoft 



momentum in Eq. i|31|) . there isn't a simple formula ex- 
pressing the correction in terms of the vacuum expecta- 
tion value of an operator involving Wilson lines like the 
one in Eq. (|23J. 

In a region \T — 1| ~ A 2 the higher order terms in the 
ultrasoft momentum that were neglected in Eq. (|31fl are 
important. Eq. (|31|) is appropriate for a region ST near 
T = 1 that satisfies 1 > 8T > A 2 , for example, ST ~ A. 



B. Jet Masses 

The squared jet masses Af 2 h are the squares of the 
invariant mass of all the particles in the two hemispheres 
a and 6, defined by the plane perpendicular to the thrust 
axis. Two commonly used variables are the sum of jet 
masses, Af| = (Af 2 + Af 2 )/Af§, and the heavy jet mass 



Ml = 



max(M 2 ,M 2 )/Ml 



The jet masses are Af 2 = 
(Pq + k { u a) ) 2 and Af 2 = (p q + k { u ] ) 2 . More explicitly, 



Mi 



Ml 



ipt+k^ + ){p- + k^-)-(p q± 



(32) 



Recall that p g is aligned along n so that p q ± = 0. Also, 
the delta function in Eq. ((PI)) sets p q ± — and pq = 
p~q = Mz- Then, working to linear order in the ultrasoft 



(29) momenta, Af 2 



M z k ( u a)+ and Af 2 



M z k u b) ~ 



dr 

dAf 2 



dr 

dMjj 



.(0) 
2-jet 



.(0) 
2-jet 



S(Ml)-S'(Ml) 
8{M 2 H ) - 5'{M%) 



M S y 

M z 

qMh 



(33) 



where 



Of-) = (max(4 a)+ ,4 b) 



(34) 



Note that in the kinematic region where expanding to 
linear order in ultrasoft and residual momentum is ap- 
propriate, the nonperturbative corrections to the Af J and 
1 — T distributions are given by the same nonperturba- 
tive matrix element. The nonperturbative corrections to 
the AfJ and M^ distributions are different. 

Working to higher orders in k u /Mz, the definitions 
of thrust in Eq. I|3U|I and of jet masses in Eq. I|32|) be- 
come different beyond linear order. However, the cor- 
rections to event shape distributions at higher orders in 
AQCD/(Afze) come not from expanding the argument of 
the delta functions used to define these variables to higher 
orders in k u /Mz, but rather from expanding these delta 
functions as power series in the ultrasoft momentum, as 
in Eq. I125[| for the jet energy. So even at higher orders, 
the enhanced nonperturbative corrections, i.e. of order 
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Pq 


- t 



FIG. 1: Determination of the thrust axis. To the order 
we are working, the quark and antiquark have momenta 
p g | = \pq\ = Mz/2. The antiquark then makes an angle 
Oq — 2 |kgx| /Mz with the z-axis, and the thrust axis t makes 
an angle 6 t = |I%±| /Mz with both the quark and antiquark. 



[Aqcd / '(Mze)] n , n > 1, come from the leading-order cor- 
rection to the argument of the delta function, which are 
the same for thrust and jet mass sum. So the enhanced 
nonperturbative corrections to thrust and jet mass sum 
are related to all orders in Aqcd /(Mze). 



C. Jet Broadening 

Jet broadening variables B ao are defined by 

xtl, (35) 



i£a,b 



where the hemispheres a and b are defined as before, and 
t is the thrust axis. The jet broadening variables at order 
k u /Mz require knowing the thrust axis to order k u /Mz- 
The thrust axis t maximizes t • P» • 

The angle between and the z-axis is given by 



q Mz/2' 
and the thrust axis t can be written as 
t = (0, — sin 6t, cos 6t). 

By symmetry, 

M z ' 



(36) 



(37) 



(38) 



which is half the size of 9q (see Fig. 

Now calculate |p* x t| for each particle. To linear order 
in k u /Mz we find for the quark, 



I il Mz ■ a 
|P<? x t| = — sm6 t 



and for the antiquark, 

Ipq x t| = |kqj_| cos6» t 



sin t> t 



(39) 



(40) 



For each ultrasoft particle a, the cross product k Q x t is 
given by the determinant 



x 

UX 





y 



z 

l,z 

sin 0t cos 9 t 



(41) 



Since sm9 t is already of order k u /Mz the cross product 
is, to linear order in k u /Mz, 



k a xt = (*s.-*S.o), 



(42) 



so |k Q x t| = |k Q j_|< Thus, to the order we are working, 
the jet broadening variables are given to this order by: 



B a = - 



E 



k a _i_| 



(43) 



where the sum on a is over the ultrasoft particles in hemi- 
sphere a or b. 

One conventionally defines two other broadening vari- 
ables as 



B n 

B s 



max (B a ,Bb) 
B a + B b 



The jet broadening distribution is 



dT 
dB 



.(0) 
2-jet 



5(B) -6'(B) 



(Of) 



M? 



(44) 



(45) 



for B Q ,fc.sum,max, where (Of) = (B) is the matrix ele- 
ment of the appropriate quantity in Eqs. (|43I44|I . Non- 
perturbative effects in the jet-broadening measures are 
not related to the jet energy or thrust. 



D. C parameter 

The C parameter is defined as 

C = 3 (A1A2 + A2A3 + A3A1) . 
where \ are the eigenvalues of 

e rs = JLy\ 

M z Z - 



P[pf 



(46) 



(47) 



and r,s — 1,2,3 are the space components of the mo- 
mentum pi of the i th particle. 

The largest component of 9 rs is 9 ZZ . The quark and 
antiquark in the jets have z-momentum p z q — — p| = 
Mz/2 to the order we are working. Then, to linear order 
in k u /Mz, the eigenvalues of 9 rs are given by: 



det(0 - XI) = (1 - A) det(AT - XI) 



(48) 
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where / is the identity matrix and 



'..a-2 



M z |k 



X V kf 



X12 — X21 



^ M z \Y a \ 



(49) 



Here the sums over a are only over ultrasoft particles. 
(The contributions from the quark and antiquark to these 
components of 9 rs are suppressed by another factor of 
1/M Z , since |p,| = |p ? | = M z /2.) 

The largest eigenvalue is Ai = 1, and the other two 
eigenvalues satisfy 



A2 + A3 = I77£ 



(KV + ( *%) 



Thus, 



The C distribution is then 



dT 
dC 



2-jet 



5(C) - S'(C) 



(o?) 



(50) 



(51) 



(52) 



where (Of) = (C) defined in Eq. (|BT|l . Like jet broad- 
ening, the C parameter distribution is not local on the 
ultrasoft fields, and the nonperturbative correction is not 
related to that for any of the above distributions. 



E. Energy-Energy Correlation and Jet-Cone 
Energy Fraction 

The angular correlations of radiated energy can be 
characterized by the one-point and two-point correla- 
tions heel 



dfldfl 



dn J <y Mz 

7 = / dTj2^-8(Sl-Sl i )S{n f -Sl j ), 



(53) 

where the sum is over all particles, and includes the terms 
with i — j. They are normalized so that 



dn 



, d 2 S _ dS 

dndw ~ d?r 



(54) 



The energy-energy correlation function P(cosx) is de- 
fined by 

P(cos X ) = f dn'dtl-^^S(cos X -cos6na>),(55) 

where 6qq/ is the angle between vectors in the f2 and Q' 
directions. 

The angular energy correlations Eq. I|53|l were defined 
111 Ref. Pfl M 

for e + e annihilation, and the solid 
angle was defined with respect to the beam direction. 
For unpolarized Z decay, there is no preferred direction, 
so dS/dfi is a constant, and d 2 £/di!d£!' contains the 
same information as the energy-energy correlation func- 
tion P(cosx)- One can, however, define distributions 
analogous to Eq. (|53|l where the solid angle is measured 
with respect to the thrust axis t. The one-point function 
is called the jet cone energy fraction J. 

The energy-energy correlation and the jet cone energy 
fraction both are proportional to S functions if one con- 
siders the leading order process of Z decay into a quark- 
antiquark pair: 



P(cosx) = ^(cosx) 
1 



(56) 



= -r [5(cosx-l) + <5(cos X + l)]. 

Ultrasoft emission (in two-jet events) changes the distri- 
bution in two ways: (a) by changing the energy or (b) 
by changing the solid angle of the emitted particles. At 
order k u /Mz, the change in energy can be neglected, be- 
cause it does not shift the angles of the partons; thus 
there is no contribution proportional to S'(cosx ±1), as 
for variables such as thrust. The angle between the quark 
and antiquark is [compare Eq. : 



cos 9 a 



k 2 

Ml 



(57) 



and the angle of the quark or antiquark with respect to 
the thrust axis is [compare Eq. J3HJ]: 



k 2 

COS 9 rl 2 = 1 — — TTT 

qt 2M% 



(58) 



where kj^ is the total _L momentum of the ultrasoft par- 
ticles. The shift in angle is second order in k u /Mz, and 
so to first order, there is no enhanced contribution near 
cosx — ±1. There are nonperturbative contributions at 
second order. 



F. Classes of Observables 

The different observables we have discussed can be di- 
vided into classes, based on the extent to which their 
nonperturbative corrections are inclusive on the ultrasoft 
degrees of freedom. 

A class I observable is the jet energy distribution. The 
nonperturbative correction to the jet energy depends on 
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where k u is the total ultrasoft momentum, so the 
jet energy distribution is totally inclusive on the ultrasoft 
fields. The derivation of nonperturbative corrections to 
the two jet energy distribution is not quite on the same 
footing as the derivation of nonperturbative corrections 
to the B meson semileptonic decay rate, because of the 
additional assumption about the equivalence of sums over 
partonic and hadronic states discussed after Eq. 10. 

Class II observables are thrust and the jet masses 
M| h . The nonperturbative corrections to these vari- 
ables require the ultrasoft momentum to be broken up 

into two parts, k u = ku + fci & ' , corresponding to the con- 
tributions from ultrasoft partons in the two hemispheres. 
The hemispheres are chosen based on the jet directions, 
i.e. based on the collinear degrees of freedom. The mo- 
mentum in each hemisphere can then be defined by inte- 
grating the ultrasoft e nergy-m omentum tensor over the 
hemisphere at infinity Hi IH HI H E3 . The class II 
variables are not totally inclusive on the ultrasoft vari- 
ables, but require them to be divided globally into two 
parts. Whether our derivation of the nonperturbative 
corrections for class II observables (e.g. the relation be- 
tween jet mass and thrust distributions) is valid depends 
on the nature of hadronization in QCD. The ultrasoft 
fields end up inside final state hadrons. The final hadron 
can contain ultrasoft partons from different hemispheres, 
so the hadronic energy flow in each hemisphere does not 
have to equal the parton energy flow in each hemisphere. 
If the hadronic and partonic energy flows differ by order 
unity, the derivation of nonperturbative effects in class 
II observables is invalid. If, for a smearing region of 
size A the mixing of ultrasoft momenta between the two 
hemispheres during hadronization is an effect of order 
Aqcd/A, then its impact on class II observables is the 
same size as fc 2 . effects, which are one higher order than 
the terms we have computed. 

Class III observables are the jet broadening measures 
B a , b, sum, max and the C parameter. These depend on 
knowing the individual ultrasoft momenta of each par- 
ton. This appears to be a notion that cannot be made 
rigorous in field theory. 



have to be imposed to reproduce other proposed relations 
amongst nonperturbative parameters for event shape dis- 
tributions HI m E . 

As an example, consider the C parameter, for which 
the nonperturbative matrix element was defined as 



(O? 




(59) 



For on-shell soft gluons collinear to the antiquark or 
quark jet (i.e. in hemisphere a or 6, respectively), fc' a ' + <C 
k^l- and fc( 6 )~ < k (b)+ . This implies that 



E 



|k Q J 2 



= 2 / V \k + k-\ 
\^ \k+ + k- 



coll 



2 £ie )+ i + Ei fc 

\ a /3 

2(e ) +e ) -) 



This leads to 



(Of) coll = 6(Of). 



(60) 



(61) 



To take into account that ultrasoft gluons can also be 
radiated at a finite angle, one can impose the condition 
that the matrix elements of Of and Of are given by the 
one-gluon contribution in perturbation theory, perform- 
ing the angular integrals in the phase space at a fixed 
value of |kj_|. Under this assumption, the matrix ele- 
ment of Of is given by: 



d.9 sin 0- 



1 



1— gluon WO 



(|ki_|/sm0) sin 2 9 



(62) 



where the factors of sin 9 from the phase space, from the 
relation |kj_| = |k| sin#, and from the squared amplitude 
for one gluon emission have all canceled out to give the 
final result. For the matrix element of Of, we calculate: 



G. Model-Dependent Relations among Event 
Shape Variables 

Nonperturbative corrections to event shape distribu- 
tions have been considered extensively in the literature 
in the past. For example, in the work of Ref. [Tl| . nonper- 
turbative shape functions were derived for thrust and jet 
mass distributions. The enhanced nonperturbative cor- 
rections to these distributions are given by first moments 
of these shape functions, and the results in sections IlII Al 
and llll Bl are in agreement with Ref. [llj . 

The derivations of the enhanced non-perturbative cor- 
rections in this section have only relied on the fact that 
they arise from matrix elements of ultrasoft operators. 
It is insightful to understand what further conditions 



fc(«)+ + fcW- ) = 2( / ABsanB 1 - 

1 -gluon \ ,/q 




(0 1 )l-gluon — - (Of )l-gluon- 



Given the assumptions that have to be made to obtain 
Eq. (|64|l (or analogous relations based on higher orders 
in perturbation theory), it does not seem likely to us 
that there is a simple analytic nonperturbative relation 
between (Of) and (Of). 
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H. Comparison with the Data 

Predictions for event shape variables have been com- 
pared with experimental data in Refs. 0,0- Nonpertur- 
bative corrections have been included using the ansatz 
that their effect on distributions for shape variables is 
described by shifting the variable by cfi np /E cm in the 
perturbative formula for the distribution. Here c is a 
constant that depends on the kinematic variable, /i np is 
a universal nonperturbative parameter, and E CTCL is the 
center-of-mass energy. An analysis in perturbation the- 
ory (similar to what was done in section Till G(l provides 
simple relations between the c's for some of the event 
shape variables. We have found that, provided one is 
not in a kinematic region that is extremely close to the 
partonic endpoint (i.e. the shape function region), c for 
1 — T and M| are the same. However, we argued that c 
for other parameters like the heavy jet mass and C are 
not connected to c for thrust. Some experimental evi- 
dence for this can be found in the analysis of Ref. Q. 
For 1 — T and the jet mass sum 2 a simultaneous fit for 
a s and /z np under the assumption that c takes on its con- 
jectured values (see Fig. 9 in 0) yields values of /x np 
that are close to each other, and values of a s that are 
consistent with other extractions of the strong coupling. 
However, Ref. Q finds that fj, np for the heavy jet mass, 
C parameter, and jet broadenings are not related to /z np 
for thrust in the way that the analysis based on pertur- 
bation theory suggests, and, furthermore, a fit to these 
variables does not yield a value of a s that is consistent 
with other extractions. 



IV. CONCLUDING REMARKS 

We have studied nonperturbative effects in Z decay 
to hadrons using soft-collinear effective theory (SCET). 
The jet energy distribution for two jet events has en- 
hanced nonperturbative effects when the jet energy is 
near Mz/2. These nonperturbative effects can be ex- 
pressed in terms of the vacuum expectation value of op- 
erators involving Wilson lines. The Wilson lines arise 
from the coupling of ultrasoft gluons to collinear degrees 
of freedom in the jet. In Appendix iBl we derive the order 

2 Ref. [3 advocates the use of a modified E-scheme jet mass to 
reduce sensitivity to hadronic masses. 



a s perturbative corrections to the jet energy distribution 
and discuss the implications of perturbative and nonper- 
turbative physics on the first moment of this distribution. 

For a region of \Ej — Mz/2\ that is of size A, the lead- 
ing nonperturbative corrections to the jet energy distri- 
bution are of order Aqcd/ A when A is large compared 
to Aqcd- In this region they can be characterized by the 
vacuum expectation value of a single operator involving 
ultrasoft fields which provides a contribution to the jet 
energy spectrum that is proportional to 5' (Mz/2 — Ej), 
For multijet events, a similar analysis holds; however, an 
additional operator analogous to 0\ but involving adjoint 
Wilson lines occurs for a gluon jet 0] . 

When A ~ Aqcd, one is in the shape function re- 
gion and the functional dependence on Ej is much more 
complicated. While we focused mostly on the kinematic 
region where Mz > A > Aqcd, it was shown that in 
the shape function region A ~ Aqcd, it is not possible 
to capture the effects of nonperturbative physics by in- 
troducing a single nonperturbative parameter /i np and 
shifting Ej — > Ej — /i np in the perturbative formula for 
the jet energy distribution. 

The jet energy distribution has the special property 
that it is totally inclusive in one of the jets, and hence ex- 
pressions for nonperturbative effects can be derived using 
operator methods that are similar to those used for the 
endpoint region in inclusive semileptonic B decay. Other 
event shape variables (for example, thrust, jet mass, jet 
broadening, etc.) have nonperturbative effects that are 
enhanced in the partonic endpoint region. We discussed 
the extent to which enhanced nonperturbative effects for 
these variables can be understood using field theoretic 
methods in QCD. 
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APPENDIX A: PROPERTIES OF WILSON LINES 

In this section we derive some useful properties of the ultrasoft Wilson lines introduced in Eq. (0 . 

1. Relations between triplet and anti-triplet Wilson lines 

Consider the time- and anti-time-ordering of the Wilson lines in the shape function S defined in Eq. Ijl6jl . For 
Y n , the path ordering is the same as time ordering and so = Y n . Consider writing Y n as the product of N 



infinitesimal integrals over path segments of length ds, 



Y, 



bjv-2 ( igAids\ b 
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(Al) 



with the subscripts on the ultrasoft gauge fields denoting their space-time location along the path of integration. 
Taking its adjoint 



yt b = fQ-igAids'j bi 



-igA N -!ds\ 



'bjv- 



bjv- 



Time ordering this expression, 
T [YL "] ■ 



^igANds\ 



'bjv-i 



bjv- 



bw_2 



l ^g— igAidsj bi 



= Y " 

1 n at 



where the overline denotes an anti-triplet Wilson line. 3 Similarly, 



T [Y na ] — Y n a , T \Y^ a ] — Y^ a 



From these results, Eq. I|18|) follows. 



1. Oi in terms of the gluon field strength 



(A2) 



(A3) 



(A4) 



We can express the operator 0\ in terms of the gluon field strength tensor as written in Eq. 124f) . It is convenient 
for this purpose to generalize the expression for the ultrasoft Wilson line to 



Y n (z; b, a) = Pexp 



ig I ds n ■ A[z + ns) 



(A5) 



so that with a = and b = oo we recover the standard Wilson line used above, Y n (z) = Y n (z; 0, oc). Differentiating 
along the n direction, 



n ■ d Yn(z) = ig I ds Yj%(z; oo, s) [n ■ d z n ■ A] (z + ns)lft (z; s, 0) 
Jo 

/•OO 

= ig / ds Yn(z\ ix, s) [n ■ d z n ■ A — n ■ d z n ■ A + n ■ d z n ■ A] (z + ns)Yn(z; s,0) 



ig / ds Y n (z; oo, s) [n ■ d z n ■ A — n ■ d z n ■ A] (z + ns)Yn(z; s, 0) 
o 

d(n-A)~ 



+ig / ds In (2;; 00, s) 
Jo 



ds 



(z + ns)Yn(z; s, 0). 



(A6) 



Using the chain rule, 



ds — [Y ri (z; 00, s) [n ■ A] (z + ns)F H (z; s, 0)] 
ds 

d , 
ds — y s (z;oo, s) 



+ / ds 1ft (2; 00, s) [n ■ A] (z + ns) 



n ■ A] (z + ns)Yn(z; s, 0) + / ds Y n {z\ 00, s) 

Jo 

^1^5,0) 



d(n • A) 



ds 



(z + ns)Y n (z; s,0) 



/ ds 1ft (z; 00, s)[n ■ A(z + ns), n ■ A(z + ns)]Yn(z; s, 0) 



+ / ds Yn(z; 00, s) 
Jo 



d(n ■ A) 



(z + ns)Yn(z; s, 0). 



(A7) 



3 Recall that the generators in the 3 representation are minus the 
transpose of those in the 3. 
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Using the above equation to eliminate the last term in Eq. ilA6(l yields 



n ■ DY n (z) = ig / dsY n (z; oo, s)n /i n 1/ G lJLV {z + ns)Y n (z; s, 0) 
Jo 

where 

n ■ DYn(z) = n ■ dY ¥l (z) — ign ■ A(oo)Y n (z) + igY n (z)n ■ A(z) 
and the gluon field strength tensor is defined by, 

G pi , = d^A u - d v A^ - ig[An,A v \. 

Hence 

1 1 r C°° 

Oi = -Tr[yJ(m -D)Y fl } = -Ti- ig / ds Y^(z; s, Q)n^n v G „ v Y n {z; s, 0) 
* ^ [ Jo 

which is Eq. l|21fl. 



(A8) 

(A9) 
(A10) 

(All) 



APPENDIX B: PERTURB ATIVE CORRECTIONS 

TO dY/dE.j 

Neglecting order a s corrections to the nonperturba- 
tive effects proportional to 0\ 1 perturbative corrections 
to dlYjct/d-Ej in Eq. © can be calculated in full QCD 
using standard methods. In this section we first review 
the computation of perturbative 0{a s ) corrections to the 
total two-jet rate r 2 _j e t and then compute the jet en- 
ergy distribution dIVjet/d-Ej at order a s . We work in 
d = 4 — e dimensions to regulate infrared, collincar and 
ultraviolet divergences that occur in contributions to the 
differential decay rate. The jets are defined using the 
Sterman- Weinberg criteria which involve an energy cut 
f3Mz and a cone half-angle 5. Corrections suppressed by 
a s p and a s S are neglected. 



1. Two Jet Decay Rate 

Using the Sterman- Weinberg definition of jets, there 
are three contributions to the two-jet rate at 0{a s ): 



(a) One quark and one antiquark each creating a jet; 



(b) One quark and one antiquark each creating a jet, 
plus a gluon with energy E g < (3Mz; 



(c) One quark and one antiquark each creating a jet, 
plus a gluon with energy E g > (3Mz inside one of 
the jets (within an angle 5 of the quark or anti- 
quark) . 

Contribution (a) is simply the rate T(Z — * qq). The 
tree and virtual gluon graphs give the amplitude: 

u s Cf 

where the color index a is summed over values a = 
1, . . . , Nc, and Cf is the Casimir of the fundamental rep- 
resentation. Explicit computation of the one-loop vertex 
correction gives, 



Mz^ q q = e^{pz)u a {p q )Y IJ 'v a {p q ) 1 + 



X ,(B1) 



4 3 2 , 
X = — r - - + ~ In 
e e 



-2Pg ■ Pg 
H 2 



4+ l2-2 ln 



-2Pq 'Pg 

M 2 



-2Pq-Pq 



Integrating the square of the amplitude over the d dimensional two body phase space gives: 



N c 
32tt 2 

x 1 



(9v + 9a) 
a s CF 



In- 



Ml 



3-e 



-4 



12 



1, 2 Ml 
2 fi 2 



3, M 2 - 
2 ln ^ 



(B2) 



(B3) 



where Qd is the total solid angle in d dimensions. The 1/e poles will cancel out against divergences from the real 
gluon emission graphs. We do not need to expand the bracketed prefactor in Eq. (|B3|) in powers of e because the 
identical factor will appear in the real gluon graphs. 
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Contributions (b) and (c) come from integrating the square of the amplitude for real gluon emission, Z — > qqg, over 
the three-body phase space in d dimensions. We find for the terms that do not vanish as (3 and 5 go to zero, 



and 



-.(c) 
Z^qqg 



Z^qqg 



g^M^NcCf 



g 2 s Mz~ e N c C F 
256?r 5 



JJ_ 

M z 



n 2 -e^3~e(gv +9a, 



^r(-S)y« 

T{-e) P ' 



(B4) 



256tt 5 



<-> 2, (4)' 



n 2 _ e n 



3-e 



3-e 



(25)" 



-(2/3)- 



/ 3e 13e 2 \r(-e) 2 
2(1 + — + 



8 /T(-2e) 



(B5) 



Adding together contributions (b) and (c) , expanding in powers of e and converting to the MS scheme yields for the 
total rate for Z — > qqg in the two-jet region 



Z^qqg 



N c C F a s 2 2 
-(9v+9a) 



32tt 3 

4 3 



Af^ £ (47r) £ ? -fi 



2. Mj 

2 



3 - e 

3, Mf I, 2 Ml 
-ln^ + -ln 2 
2 /i^ 2 /i 



(B6) 



+ - - - In - ~ In + ~ ln z - 41n2/31n<5 - 31n<5 + ^ 



13 llTT 2 

2~ ~ 12 



Finally, we add together the rates Y z-* q q and r 'z^qqg from Eqs. I|B3|) and l|B6|) . The e-dependent prefactors in brackets 
in the two equations are identical, as promised. The 1/e-poles in the remainder of the expressions cancel out exactly 
(as do all the logarithms of Mz/fJ>), so we can set e = in the remaining finite parts, leaving 



NcMz, 2 , 2N 

2-jct = — s — (g v +9 a) 
iZn 



1 + asCFr*_^_ 3ln§ _ Unmn§ 

7r V 2 3 



(B7) 



which agrees with Sterman and Weinberg's original result |2^ . 



2. Differential Decay Rate dr 2 -j e t /dEj 



We now turn our attention to the differential decay rate dlVjet/d-Ej. The contribution of Tz—>qq to this rate is 
simply 



dT Z^qq _ r M Z 



(B8) 



where Tz^qq is the total rate for Z — > qq calculated to 0(a s ), which is given by Eq. I|B3(1 . 

For the contribution of real gluon emission processes, we write the three-body phase space for this rate: 



dr 



2-jct 



1 



1 



fl d ^d-idE 1 E^- 4 dE 2 E^-' i d cos 9 sin a " 4 6 



d-4 



dEj 



16M Z (2tt) 



2d-3 



(B9) 



xS 



M 2 z -2M z {E l +E 2 
2 Ei E2 



1 — cos t 



S(Ej 



where the S(Ej — ■ ■ •) defines Ej according to which par- 
tons actually go inside the jet. It is useful to split up the 
phase space slightly differently than for the case of the 
total rate: 

(a) Gluon with energy E g > (3Mz inside unobserved 
jet; 

(b) Gluon with any energy inside observed jet; 

(c) Gluon with energy E g < f3Mz outside observed jet. 



These three regions exhaust the possible gluon energies 
and locations with respect to the jets. It is convenient to 
introduce the variable 

M 7 

ej = ^-Ej (BIO) 

and focus on a region of ej near the origin with size of 
order (3M Z - 

For case (a), where a gluon with E g > f3Mz is inside 
the unobserved jet, take E\ = E g , E% = Eg, so 9 is the 
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angle between the gluon and antiquark, and Ej = E q . 
Integrating over 9 and Eg using the delta functions leaves 
an integral over E g running between the limits 



and restricts e.j to lie between 



S 2 /3M Z < ej < 



M Z S 2 



(Bll) 



(B12) 



Similarly, for case (b), where a gluon with any energy 
lies inside the observed jet, E\ — E q , E2 = E g , and 
Ej = Eg + E q . Integrate over 9 and E q using the delta 
functions. Then the limits of the E g integral are 



E, 



M, 



i± Wi 



4 V" V ' M Z S 2 
and ej is restricted to the region 



M Z S 2 



< ej < 0. 



(B13) 



(B14) 



Before we proceed to case (c), note that the physical 
observable we actually want to calculate is the smeared 
distribution 



dr 



dE, 



= J dE jWA (Ej)-^-, (B15) 



where u>a is a smooth function which smears the dif- 
ferential rate over a region of jet energy whose size is 
of order (3Mz- But the contributions to the rate from 
cases (a) and (b) have support only over a region of size 
5 2 Mz <C (iMz near ej = 0. Consider smearing dT/dEj, 
or equivalently dr/dej, over a region near Ej — Mz/2 
(ej = 0) of size of order j3M z . Then w(0) ~ l//3M z , 
w'(0) ~ 1/((3M Z ) 2 , etc. Expanding, 

/dr f dr 
dejw(ej)— = / de,j[w(0) + w'(0)ej + ■■■} — . 
dej J dej 

(B16) 

Since w 1 (0)/w(0) ~ l//3Mz, and, for the contributions in 
cases (a) and (b), ej ~ S 2 Mz in the region where dT/dej 
is nonzero, the second term is suppressed by a power of 
S 2 1 (3 -C 1. Thus only the first term is relevant. 4 Keeping 
only the first term amounts to replacing the full dT/dej 
by 



dr 

dej 



S(ej) / de'„ 



dr 
r dZ 



(B17) 



However, integrating the contributions of (a) and (b) to 
dr /dej over all allowed values of ej simply gives their 



This argument assumes that the integral J dej ejdr/dej is fi- 
nite, which can easily be shown. 



contribution to the total Sterman- Weinberg jet rate, that 
is, they will build up part of the term T z-> q q g 5{Ej — 
Mz/2) in dlVjet/d-Ej. Since we have already calculated 
the total rate, we need not analyze cases (a) and (b) any 
further, as long as we can get the remaining contribution 
to the total rate from case (c). 

In case (c) we have a gluon with E g < f3M anywhere 
outside the observed jet. Here E\ = E g , E2 = E g , and 
Ej = E q . Writing out the formula for the rate explicitly, 



dr (c) . 1 1 

dEj - 16M Z (2tt)2<*-3 d - 1 d ~ 1 
x9{ej)9[pMz{l-5 2 )~ej] 



(B18) 



ej(l+5 2 ) 



dE g E d g - i E d g-U^9\M^gg g \ 2 . 



The part of the amplitude that gives a contribution that 
survives as (3 — > is 



M 



(c) 

Z^qqg 



2 -AN r -2.,e rf - 2 ^2 , MzPq-Pq 



(k-p q )(k-Pq) 

(B19) 

Substituting Eq. (|B19I) into the phase space 



(B20) 



dT^) = M z g 2 N c C F / fi \ e 2t d-2 
dEj 256tt 5 \M Z ) 7T> d-1 



x9[f3M z (l - S 2 ) - ej] In 



l3Mz 



ejS 2 



The factor, (l/ej)\n[((3Mz — ej)/(ejS 2 )], is singular as 
ej — > 0, and must be rewritten in terms of an integrable 
quantity. Use the "plus distribution" : 

rPMz rl3M z 

dej f(ej)+ g(ej) = / dej f{ej)[g{ej) - g(0)], 

IQ JO 

(B21) 

where / diverges at ej = and g is a test function finite 
at ej = 0. To replace / by /+, we would write 



0M Z rPM z 

dejf(ej)g(ej) = / dej f(ej)+ g{ej) (B22) 
(i Jo 

rPMz 

+5(0) / dej/(ej). 



The second term amounts to replacing 

f(ej)^S(ej) [ dejf(e'j). 



(B23) 



But making this replacement in Eq. i|B20f> means writ- 
ing a delta function 5(Ej — Mz/2) and integrating the 
differential rate over all allowed values of ej, which 
again just gives its contribution to the total Sterman- 
Weinberg jet rate. Together with the contributions 
from (a) and (b) this gives the one loop contribution 
to S(Ej — Mz /2)T2-jct- Only the plus function piece 
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gives a deviation of the jet energy distribution away from 
Ej = Mz/2. The final result for the differential rate to 
0{a s ) is: 



dr 



2-jct 



dEj 



+ Mza ° 2 N f F { 9 l + 9 \) 

x9(ej)e(/3Mz - ej) 



(B24) 



llnf 



(3M 2 



ej 



ej \ ejS 2 
where the total rate r 2 _j e t is given by Eq. (|B7(I . 

3. First Moment of the Jet Energy Distribution 



0.04 



0.03 - 



0.02 



0.01 



0.2 



0,-1 



0.6 

r 



0.B 



As an application of the above result consider the first 
moment of the jet energy distribution, defined by 



Mi(f) 



M z 



dEj 

— fWz 



1 



dr 2 -jot 



r 



2-jet 



dE, 



Ej_ 
M z 



(B25) 

Using the expression in Eq. 121J) for the nonperturbative 
correction and in Eq. (|B24|) for the order a s perturbative 
correction to the jet energy distribution gives 



Mi(f) 



o: s Cf(3 



flog (jp) -(l-/)log(l-/) 



Mz ' 



for / < 1 and 



„, ,ns ot s C F f3 ( 1 
M l( /) = — — log - 



)|Oi|0) 



(B26) 



(B27) 



for / > 1. Note that the order a s contribution to M 1 (f) 
is independent of / for / > 1. This occurs because the 
perturbative correction vanishes for Ej < Mz/2 — (3Mz- 
In Fig.[21we plot Mi(/), for / < 1. For this figure the 
value of the energy cut is (3 = 0.15 and the cone half- 
angle is 5 = 15° and the vacuum expectation value of 
0\ is set equal to 500 MeV. We evaluate a s at the scale 
(3Mz and find with these parameters that the order a s 
corrections reduce the two jet rate by about 16% from its 
tree level value lending support to the validity of pertur- 
bation theory for the values of the cone angle and energy 
cut used in Fig. 



FIG. 2: Plot of the function Mi(/). The black solid curve 
shows the perturbative contributions only, while the red 
dashed line represents the moment including the nonpertur- 
bative contribution. The figure corresponds to j3 = 0.15, 
S = 7r/12, (0|Oi|0) = 0.5 GeV, and we have evaluated the 
strong coupling constant at the scale fi = (3Mz- 



The full theory amplitude for Z — > qq is reproduced in 
SCET by the Wilson coefficient in the current matching: 



JQCD = [^Wn)C(n,p q -p q )T^[W^ n ] 



(B28) 



where there is an implicit sum over label momenta and 
the matching coefficient C(/i, p q -pq) can be read off 5 from 
Eqs. IE) and (|B2|) : 



C(n,p q -Pq) 



a s C F 



2n 



12 



1,2 {- 2 Pq-Pq 



2 



-2p q - Pq 



(B29) 



and the UV renormalization factor for the current in the 
effective theory is 



2tt 



4 

72 



e V M 



(B30) 



Note that both the renormalization factor and the match- 
ing coefficient depend on the label momenta for the quark 



4. Perturbative Corrections in the Effective 
Theory 

Although we have used full QCD to calculate the jet 
energy distribution it is possible to do the computation 
in the effective theory. Here we briefly discuss how that 
computation would proceed. 



5 The matching coefficient is just given by the finite part of the 
full theory matrix element (gg|j M |0) because the full theory cur- 
rent has no anomalous dimension, so the 1/e poles are pure IR 
divergences, which must cancel out in the matching condition. 
The loop graphs in the effective theory contributing to this ma- 
trix element are zero in dimensional regularization, so the finite 
part of the matching coefficient is just the finite part of the QCD 
matrix element, given by Eqs. |EJ and iB2l . while the infinite 
parts become the U V counterterm in the effective theory l.iii l37l . 
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and antiquark. For outgoing particles, the collinear Wil- 
son lines are defined as 



W n (z)=Pexp 



ds n ■ A n (ns + z) 



(B31) 



and one must include collinear gluons produced by a Wil- 
son line in real gluon emission to get the correct Z — > qqg 
amplitude. We find that the perturbative expressions for 



the two jet rate presented in the previous sections are 
reproduced by the effective theory if we call any par- 
ticles inside the observed quark jet n-collinear particles 
and all other particles n-collinear. In the effective theory, 
ultrasoft gluons in the final state contribute zero in per- 
turbation theory and appear only in the nonperturbative 
shape function. A similar result holds for deep inelastic 
scattering 38]. 
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